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CN ■ Abstract 



The extended Lomon-Gari-Krumpclmann model of nucleon electromagnetic form factors, 
\ which embodies p, p' , u, lo' and (j> vector meson contributions and the perturbative QCD 

high momentum transfer behavior has been extended to the time-like region. Brcit-Wigncr 
Cl^ formulae with momentum-dependent widths have been considered for broad resonances in 

order to have a parametrization for the electromagnetic form factors that fulfills, in the 
time-like region, constraints from causality, analyticity, and unitarity. 

This analytic extension of the Lomon-Gari-Krumpclmann model has been used to perform 
a unified fit to all the nucleon electromagnetic form factor data, in the space-like and time-like 
region (where form factor values are extracted from e + e~ O NN cross sections data), 
^v^j \ The knowledge of the complete analytic structure of form factors enables predictions at 

7-H ■ extended momentum transfer, and also of time-like observables such as the ratio between 

^sO | electric and magnetic form factors and their relative phase. 



^ ! 1 Introduction 



X 

5-H 



Nucleon electromagnetic form factors (EMFF's) describe modifications of the pointlike photon- 
nucleon vertex due to the structure of nucleons. Because the virtual photon interacts with single 
elementary charges, the quarks, it is a powerful probe for the internal structure of composite 
particles. Moreover, as the electromagnetic interaction is precisely calculable in QED, the dy- 
namical content of each vertex can be compared with the data. 

The study of EMFF's is an essential step towards a deep understanding of the low-energy QCD 
dynamics. Nevertheless, even in case of nucleons, the available data are still incomplete. 
The experimental situation is twofold: 

• in the space-like region many data sets are available for elastic electron scattering from 
nucleons (N), both protons (p) and neutrons (n). Recently, the development of new 
polarization techniques (see e.g. Ref. pQ) provides an important improvement to the 
accuracy, giving a better capability of disentangling electric and magnetic EMFF's than 
the unpolarized differential cross sections alone. 

• In the time-like region there are few measurements, mainly of the total cross section (in a 
restricted angular range) of e + e~ <-> NN, one set for neutrons and nine sets for protons, 
one of which includes a produced photon. Only two attempts, with incompatible results, 
have been made to separate the electric and magnetic EMFF's in the time-like region. 
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Many models and interpretations for the nucleon EMFF's have been proposed. Such a wide 
variety of descriptions reflects the difficulty of connecting the phenomenological properties of 
nucleons, parametrized by the EMFF's, to the underlying theory which is QCD in the non- 
perturbative (low-energy) regime. 

The analyticity requirement, which connects descriptions in both space (q 2 < 0) and time-like 
(q 2 > 0) regions, drastically reduces the range of models to be considered. In particular, the 
more successful ones in the space-like region are the Vector-Meson-Dominance (VMD) based 
models [2j[3] (see, e.g., Ref. [I] for a review on VMD models) that, in addition, because of their 
analytic form, have the property of being easily extendable to the whole g 2 -domain: space-like, 
time-like and asymptotic regions. 

In this paper we propose an analytic continuation to the time-like region of the last version 
of the Lomon model for the space-like nucleon EMFF's [5]. This model has been developed 
by improving the original idea, due to Iachello, Jackson and Lande [2] and further developed 
by Gari and Krumpelmann [3], who gave a description of nucleon EMFF's which incorporates: 
VMD at low momentum transfer and asymptotic freedom in the perturbative QCD (pQCD) 
regime. As we will see in Sec. EJ in this model EMFF's are described by two kinds of functions: 
vector meson propagators, dominant at low-g 2 and hadronic form factors (FF's) at high-g 2 . The 
analytic extension of the model only modifies the propagator part and consists in defining more 
accurate expressions for propagators that account for finite-width effects and give the expected 
resonance singularities in the g 2 -complex plane. 

2 Nucleon electromagnetic form factors 

The elastic scattering of an electron by a nucleon e~N — > e~N is represented, in Born ap- 
proximation, by the diagram of Fig. Q] in the vertical direction. In this kinematic region the 
4-momentum of the virtual photon is space-like: q 2 = —2^1^2(1 — cos# e ) < 0, 0J\^2) is the 
energy of the incoming (outgoing) electron and 8 e is the scattering angle. 




Annihilation 



Figure 1: One-photon exchange Feynman diagram for 
scattering e~N — > e~N and annihilation e + e~ — > NN. 

The annihilation e + e~ — > NN or NN — > e + e~ is represented by the same diagram of Fig. Q] but 
in the horizontal direction, in this case the 4-momentum q is time-like: q 2 = (2co) 2 > 0, where 
lo = uj\ = 0J2 is the common value of the lepton energy in the e + e~ center of mass frame. 
The Feynman amplitude for the elastic scattering is 

M = \[eu(k 2 )^u(h)] [eU( P2 )T fl (p 1 ,p 2 )U( Pl )] , 

where the 4- momenta follow the labelling of Fig. [H u and U are the electron and nucleon spinors, 
and is a non-constant matrix which describes the nucleon vertex. Using gauge and Lorentz 
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invariance the most general form of such a matrix is [6] 



(i) 



where Mjy is the nucleoli mass (N = n, p), Ff(q 2 ) and F? (q 2 ) are the so-called Dirac and 
Pauli EMFF's, they are Lorentz scalar functions of q 2 and describe the non-helicity-flip and the 
helicity-flip part of the hadronic current respectively. Normalizations at q 2 = follow from total 
charge and static magnetic moment conservation and are 



F l N (0) = Q N , FJ?(0) 



N i 



where Qn is the electric charge (in units of e) and hn the anomalous magnetic moment (in units 
of the Bohr magneton of the nucleon N. 

In the Breit frame, i.e. when the transferred 4-momentum q is purely space-like, q = (0,q), the 
hadronic current takes the standard form of an electromagnetic 4-current, where the time and 
the space component are Fourier transformations of a charge and a current density respectively: 



Pq 



AM 2 



F 



N 



J q = eU(p 2 )iU(p 1 )[F 1 N + F 2 N ] 



We can define another pair of EMFF's through the combinations 

„2 



AM 



2 F 2 

N 



N 



(2) 



Ff + F 2 N 



these are the Sachs electric and magnetic EMFF's [7], that, in the Breit frame, correspond to the 
Fourier transformations of the charge and magnetic moment spatial distributions of the nucleon. 
The normalizations, which reflect this interpretation, are 



G»(0) = Q N , GZ(0)= m , 

where = Qn + K N is the nucleon magnetic moment. Moreover, Sachs EMFF's are equal to 
each other at the time-like production threshold q 2 = AM^, i.e.: 



G%(AM 2 N ) = GU±M%). 
Finally, we can consider the isospin decomposition for the Dirac and Pauli EMFF's 

- FP) , i = 1, 2 , 



(3) 



F 1S and F lv are the isoscalar and isovector components. 



3 The Model 

The model presented here is based on simpler versions designed for the space-like EMFF's of 
Iachello, Jackson and Lande [2] and of Gari and Krumpelmann idea [3], which describes nucleon 
EMFF's by means of a mixture of VMD, for the electromagnetic low-energy part, and strong 
vertex FF's for the asymptotic behavior of super-convergent or pQCD. The Lomon version [5], 
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which fits well all the space-like data now available included two more well identified vector 
mesons and an analytic correction to the form of the p meson propagator suitable for describing 
the effect of its decay width in the space-like region fitted to a dispersive analysis by Mergell, 
Meissner, and Drechsel [8]. 

This model describes the isospin components, eq. ([3]), in order to separate different species of 
vector meson contributions. For the isovector part the Lomon model used the p and /?(1450) or 
p' contribution, while for the isoscalar the lj, w(1420) or lj' and (ft were considered. In detail 
these are the expressions: 

F?{q 2 ) = [BW^ B (q 2 ) + BWgtf)]Ff (q 2 ) + [l - BW^ D (0) - BWg '(0)]jf V) 
F?{q 2 ) = k p [BW^ B (q 2 ) + Bwg V)] W) + [l " BW^(0) - BW$ '(0)] if \q 2 ) 

Ff(q 2 ) = [BW£(q 2 ) + BW^' (q 2 )]F? (q 2 ) + BW$tf) Ff(q 2 )+ (4) 
[l-BWtf(0) - BW$'(0)]FV(q 2 ) 

Ff{q 2 )=[K i0 BW^{q 2 ) + k u , BW^'(q 2 )]F^ \q 2 ) + ^ BW$ \q 2 ) 1$ \q 2 ) + 
[ks - k u BW$(0) - k u , BW$'(0)]F 2 D (q 2 ) , 



where: 

• BW$(q ) is the propagator of the intermediate vector meson a in pole approximation 
BW Q a (q 2 ) = 9 f , a = p',u;,u/,0, 

M 2 g a / f a are the couplings to the virtual photon and the nucleons; 



^^MMd(? 2 ) are dispersion-integral analytic approximations for the p meson contribution 
in the space-like region [8] 

1)P 2 _ 1.0317 + 0.0875 (1 - g 2 /0.3176)- 2 
^mmdWJ- 2(1-^/0.5496) 

r%P ,J2\ _ 5 - 7824 + °- 3907 i 1 ~ gVO-1422)- 1 



2kJ1 - g 2 /0.5362) 



the last term in each expression of eq. (jlj) dominates the asymptotic QCD behavior and 
also normalizes the EMFF's at q 2 = to the charges and anomalous magnetic moments 
of the nucleons; 

the functions F®(q 2 ), a = p,lj,4> and i = 1,2, are meson-nucleon FF's which describe the 
vertices aNN, where a virtual vector meson a couples with two on-shell nucleons. Noting 
that the same meson-nucleon FF's are used for p' and lj' as for p and lj, we have 

/ 2 \ 3/2 

*?V) = /i(? 2 )(^Afj » (6) 



3/2 

Fi{q 2 ) = f2(q 2 )' 



Af q 2 -4\ 
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where Ai and A2 are free parameters that represent cut-offs for the general high energy 
behavior and the helicity-flip respectively, and 



- 2 _ 2 m[(Aj,-a/A 2 QCD ] 

" " ln(A^/A^ CD ) ' W 

where Ad is another free cut-off which controls the asymptotic behavior of the quark- 
nucleon vertex, the extra factor in Ff(q 2 ) imposes the Zweig rule; 

the functions FP(q 2 ) can be interpreted as quark-nucleon FF's that parametrize the direct 
coupling of the virtual photon to the valence quarks of the nucleons, 

F?(* 2 ) = JK.(tA^)\ -1,2, (8) 



Al-cp \A 2 -q 2 

q 2 is defined as in eq. ([7]); 

• finally, n a is the ratio of tensor to vector coupling at q 2 = in the aNN matrix element, 
while the isospin anomalous magnetic moments are 

K, s — Kp -\- K n , K v — Kp K n . 

The space-like asymptotic behavior (q 2 — > —00) for the Dirac and Pauli EMFF's of eq. (jJJ) is 
driven by the Ff 2 (Q 2 ) contribution, given in eq. ([8]). In particular we get 



Fr\q 2 



1 



q 2 ln(-q 2 /A 2 QCD ) 

Fi v ' is (q 2 ) 
- q Hn(-q 2 /A 2 QCD ) 



as required by the pQCD [9]. 

In principle this model can be extended also to the time- like region, positive q 2 , to describe 
data on cross sections for the annihilation processes: e + e~ f-> NN. However, a simple analytic 
continuation of the expressions given in eq. (Jl|) involves important issues mainly concerning the 
analytic structure of the vector meson components of the EMFF's that, in the time-like region, 
are complex functions of q 2 . The hadronic FF's of eqs. (J6j) and ([8|) may also have real poles as a 
function of q 2 . In fact as defined above Ff has a real pole at q 2 = A 2 . In the other denominators 
of eqs. ([6]) and (|8|), as in F^ u and Ff, q 2 is replaced by q 2 . The latter as a function of q 2 has a 
maximum in its real range < q 2 < A 2 D , which, for reasonable values of Ad and Aqcd> may be 
smaller than A 2 , A| and A 2 D . Therefore all the hadronic FF's real poles may be avoided by also 
replacing q 2 by q 2 in the factors of Ff . This does not effect the asymptotic behavior required 
by the Zweig rule and will be adopted in the model used here. The results in Sec. [6] show that 
with this modification real poles can be avoided in every case examined, although in half the 
cases mild constraints on Ai or Aqcd are needed which affect the quality of the fit negligibly. 
A detailed treatment of the possibility of extending the model from the space-like to the time-like 
region, will be given in Sec. [5j 
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4 Analyticity of Breit-Wigner formulae 



The standard relativistic Breit-Wigner (BW) formula for an unstable particle of mass M and 
energy independent width T is 

BW(s] 



M 2 -s-iTM' 

it has a very simple analytic structure, only one complex pole and no discontinuity cut in its 
domain. Once this formula is improved to include energy dependent widths one immediately 
face problems concerning the analyticity. 

We consider explicitly the case of the p resonance in its dominant decay channel 7r + 7r _ . A 
realistic way to formulate an energy dependent width is to extend the p mass off-shell, making 
the substitution M 2 = s, in the first order decay rate 

\ g P^ M 2 p - 4Ag 

where g^n is the coupling constant and, M p and M n are the p and pion mass respectively. Such 
a decay rate has been obtained by considering, for the vertex / 07r + 7r~, the pointlike amplitude 

where is the polarization vector of the vector meson p, and p± the 4- momentum of tt^. 
Finally, assuming the tt + tt~ as the only decay channel and using eq. © for the corresponding 
rate, the energy dependent width can be defined as 

T"(c)-T" M '( S ~ S ° iS - S0)l 7 ~ rgM ' (10) 

si) ~ ° * \Mi p -sJ -M p s ' 7 '- W _ fl0 )i' ( } 

where the subscript "s" indicates the factor 1 /s appearing in the width definition, Tq is the total 
width of the p, and sq = AM 2 . It follows that the BW formula becomes 



BW s (s) 



3 • 



s{M 2 -s)-i lp (s-s ) 5 



In this form the BW has the "required" |10j discontinuity cut (so,oo) and maintains a complex 
pole s p 



I p 2 + ifp M p ~ M p 2 + iT p M p = s° p 



Due to the more complex analytic structure the new pole position s p turns out to be slightly 
shifted with respect to the original position S®. Moreover, these are not the only complications 
introduced by using T p (s) instead of Tq, the power 3/2 in the denominator and the factor 1/s, see 
eq. ()10p . generate also additional physical poles which, in agreement with dispersion relations, 
must be subtracted, as discussed below. 

4.1 Regularization of Breit-Wigner formulae 

We consider the general case where there is a number iV of poles lying in the physical Riemann 
sheet. We may rewrite the BW by separating the singular and regular behaviors as 

BW(s) = _ 



Yl» =1 ( s - Zj )[M 2 -s-i 7 (s-s ) 
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where Pn(s) is a suitable N degree polynomial, j3 is a non-integer real number which defines 
the discontinuity cut (in the previous case we had f3 = 3/2), 7 = MT°/{M 2 - s f, and the 
Zj are the real axis (physical) poles. To avoid divergences in our formulae, we may define a 
simple regularization procedure consisting in subtracting these poles. In other words we add 
counterparts that at z = Zj behave as the opposite of the i-th. pole. In more detail, we may 
define a regularized BW as 

BW(s) = BW{s) - Y . (11) 

ti nf=i,;^ - *i) [M 2 - Zk - i 7 (z k - s )P] s - z k 

In the Appendix A we show how dispersion relations (DR's) offer a powerful tool to implement 
this procedure without the need to know where the poles are located. However in this paper we 
show that an analytic expression also contains the information. 

4.2 Two cases for T(s) 

In our model for nucleon EMFFs, widths are used only for the broader resonances: p(770), 
p(1450) and w(1420) [11] . We consider explicitly two expressions for T(s) which entail different 
analytic structures for the BW formulae. Besides the form we discussed in Sec. 01 eq. (jlpp . we 
consider also a simpler expression (closer to the non-relativistic form) , hence for a generic broad 
resonance we have 

M 2 ( s-~s \ I _ 7s (f ~ s~o) f ... T M 3 



s \M*-s j M s (M 2 -s )^ 

s o V - 71 / - a r M 



(12) 



r i( g ) = r M 2 - =^( S - S ~0)^ With: ^l 

M ~ s o J M (M 2 -I ) 



In both cases we assume that such a resonance decays predominantly into a two-body channel 

whose mass squared equals sq. The subscript "1" in the second expression of eq. (|12p indicates 

that there is no extra factor 1/s in the definition of the energy-dependent width. 

As already discussed, the BW formulae acquire a more complex structure as functions of s, as 

a consequence unwanted poles are introduced. Such poles spoil analyticity, hence they must be 

subtracted by hand or, equivalently, using the DR procedure defined in Appendix A. 

More in detail, for both BW formulae we have only one real pole, that we call s s and si 

respectively (both less than sq). The corresponding residues, that we call R s ,i, are 



R 



M 2 - 2s s + §7 s Vso - s s ' 
Ri = — : s 7= ■ 

-1 + 271VS0 - s l 

Following eq. (fTTj) . the regularized BW formulae read 

s - s S) i 

In particular, below the threshold sq, where BW's are real, we have 

s R s 



(13) 



BW s {s < Sq) 
BWi(s < s ) 



s(M 2 -s)-7 s (s -s) 3 / 2 s- 
1 Rt 



(14) 



M 2 — s — 71 (so — s) 3 / 2 s — si 
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Above so BW's become complex, real and imaginary parts are obtained as limit of BW s \{s) 
over the upper edge of the cut (s*o, oo). Since the poles s s> \ are real only the real parts have to 
be corrected as 

s 2 (M 2 - s) 



Re 



Re 



BW s {s > s ) 
BWi(s > s ) 



Rs 

s 2( M 2_ s )2 +7 2( s _g )3 s _ Sg 

M 2 — s ft 



(M2-s)2 +7 2( s _5 )3 s _ Sl 

while the imaginary parts remain unchanged 

SJs(s - S ) 3/2 



Irn 



Im 



BW s {s > s ) 
BWi(s > s ) 



s 2( M 2_ s )2 +7 |( s _ 5o )3 ' 

7i( S "5o) 3/2 



( M 2-s) 2 + 7i 2 (s-so) 3 



Resonance 


M (GeV) 


T (GeV) 


«o 


s s (GeV 2 ) 


si (GeV 2 ) 


p(770) 


0.7755 


0.1491 


4M 2 


0.005953 


-11.63 


p(1450) 


1.465 


0.400 


4M 2 


0.003969 


-29.43 


w(1420) 


1.425 


0.215 


(M n + M p ) 2 


0.06239 


-19.46 



(15) 



(16) 



Table 1: Parameters for the BW formulae of resonances: /3(770), p(1450) and a; (1420). 



The parameters of the subtracted poles for the three vector mesons are reported in Table [TJ A 
third case is discussed in Appendix B. It is not fitted to the data because its resonance structure 
is intermediate between the two above cases. 




(GeV) 




-2 

s 1/2 (GeV) 



Figure 2: Left: comparison between the two descriptions of the p peak in the space-like region. 
Right: relative differences w.r.t. BW S . 



Figures [2] and [3] show comparisons between the two descriptions in case of p in the space-like 
and time-like regions respectively. On the left of each figure is the modulus of each, on the right 
the relative difference with respect to BW S . 



S 




5 The analytic extension 

The original model, described in Sec. [3] and constructed in the space-like region, can be analyt- 
ically continued in the time-like region using the regularized BW formulae obtained in Sec. [H 
We consider then a new set of expressions for F" 2 (q 2 ) and Ff 2 (q 2 ), homologous to those of 
eq. where now we use regularized BW formulae instead of the MMD [8] p width form or 
the zero-width approximation given in eq. ([5]), and also two additional vector meson contribu- 
tions, p(1450) and cj(1420) here simply p' and u' , as in the last version of the Lomon model [5]. 
Such BWs have the expected analytic structure and reproduce in both space-like and time-like 
regions the finite-width effect of broad resonances. The narrow widths of the uj and eft have 
negligible effects, so we use these modified propagators only for broader vector mesons, namely: 
the isovectors p and p' , and the isoscalar lo' . These are the new expressions for the isospin 
components of nucleon EMFF's 

F'l^M 2 ) = [BWLe(Q 2 ) + BwLe(Q 2 )]F((Q 2 ) 

+ [1 - ML(0) - Bfi(0)]^( f 2 ) 

+ [k v - k p BW P caae (0) - V Bwi(O)] F 2 D (q 2 ) 

FlUM 2 ) = [BW^q 2 ) + BW^ c (q 2 )]Fr(q 2 )+BW^q 2 )Ff(q 2 )+ 
[1 - BW$(0) - BwtsM]F?{q 2 ) 

^casefe 2 ) = [^BW^(q 2 ) + K w/ BW W c UQ 2 )]F^(q 2 )+K (l} BW^q 2 )F^q 2 ) + 
[k s - k w BW$(0) - Wl(0) - K 4BWf(p)]F?tf) , 
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where case=s and case=l correspond to the parametrizations of the energy dependent width 
described in Sec. 14.21 Following eqs. f|14j) - (|16[) for the definition of BW(q 2 ), and including the 
coupling constants, we have 



BW (a 2 




l )-i^s{q 2 -4f/ 2 



R 



ill ("? 2 



case 



case 



with: (3 = p, p', u) 1 (parameters in Table [T]) and where: the 7^ are given in eq. (|12p and the 



residues R± in eq. (|13p . The introduction of the regularized BW's does not spoil the high 

energy behavior of the resulting nucleon EMFF's. In fact, as \q 2 
vanishes like 1/q 2 , i.e. following the same power law as the previous BW^{q 2 
and case=s, see eq. ©, indeed we have 



00, the function BW^ casc (q 2 ) 
in the case=l 



BW case (q 2 ) 



9pM 2 l _ Rt 



9f>M$ R{ 



case 



case 



It is interesting to notice that in both cases is just the subtracted pole which ensures the expected 
behavior and, in particular, the asymptotic limit of: q 2 ■ BW case (q 2 ) is proportional to and 
(1 — Rs) respectively. For the reason discussed at the end of Sec. El for the present model q is 
replaced by q in the hadronic FF's Ff of eq. ©. 



6 Results 

Nine sets of data have been considered, six of them lie in the space-like region [12] and three in 
the time-like region |13H21| . The data determine the Sachs EMFF's and their ratios. The fit 
procedure consists in defining a global x 2 as a sum °f nine contributions, one for each set. More 
in detail, we minimize the quantity 

9 

2 \ " 2 

i=l 

where the coefficients Tj weight the i th contribution, we use r, = 1 or t\ = to include or exclude 
the i th data set. The single contribution, xh ^ s defined in the usual form as 

H J ' 

where Qi(q 2 ) indicates the physical observable, function of q 2 , that has been measured and the 
set {ql,v l k , 6Vf.; Ni} represents the corresponding data; v\ is the k th value (k = 1, . . . , Ni) of the 
quantity Qi (i = 1, . . . , 9) measured at q 2 = q\, with error 5v l k . 

Table [2] reports the complete list of observables, the number of data points and the corresponding 
minimum x 2 ' s ! m the two considered cases as described in Sec. [5] for the sets of data with and 
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D 


TV- 


minimum Xi 


case = s 
With BaBar 


case = 1 
With BaBar 


case = s 
No BaBar 


case = 1 
No BaBar 


space-like 


I VI 


68 


48.7 


50.1 


54.6 


60.8 




36 


30.4 


27.6 


26.2 


35.0 


M 


65 


154.6 


154.2 


158.2 


167.0 


rvn 


14 


22.7 


23.2 


24.1 


26.0 


fi p G p E /G p M 


25 


13.9 


12.9 


10.6 


14.4 




13 


11.3 


10.7 


8.2 


8.9 


time-like 


\r<P i 

l^effl 


81 (43) 


162.5 


166.7 


62.2 


35.0 


i i 

l^effl 


5 


8.4 


6.3 


3.2 


0.3 




Total 


313(275) 


452.5 


451.7 


347.3 


347.4 



Table 2: Measured quantities, numbers of data points and x 2 contributions. The values in 
parentheses indicate the number of data points in the case "No BaBar" . 



without the BaBar data which have a final state photon. For case=s, with and without the 
BaBar data, the optimization over the full set of 13 free parameters (Table [3]) determines Ai, 
A2, Ajj and Aqcd such that the hadronic FF's have no real poles. For the case=l with BaBar 
data the full minimization implies a zero for (Aj — q 2 ) producing poles in the hadronic FF. Re- 
minimizing with the constraint Ai = 0.5 GeV, just above the 0.4744 GeV obtained without the 
constraint, removes the poles. For case=l without BaBar data it is required that the already 
fixed Aqcd = 0.15 GeV be changed to Aqcd = 0.10 GeV to avoid a zero of (A^ — q 2 ). In both 
cases the change in \ 2 is negligible. 

Data and fits, black and gray curves correspond to case=l and case=s respectively, are shown 
in Figs. [MEl In the space-like region the electric Sachs EMFF's are normalized to the dipole 
form 



G D (q 2 



0.71 GeV' 



while magnetic EMFF's are also normalized to the magnetic moment. This normalization de- 
creases the range of variation, but the curves clearly demonstrate deviations from the dipole 
form. The observable Rn is defined as the ratio Rn = G^/G^ for the nucleon N. As N stands 
for both neutron and proton there are six space-like observables. A departure from scaling is 
shown in the deviation of R p and R n from unity. 



The time- like effective FF, 



I G^ s I , is defined as 



cne^e 



Ana -2 



NN) 



1/2 



(17) 
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where a(e + e~ — > NN) is the measured total cross section and the kinematic factor at denomina- 
tor is the Born cross section for a pointlike nucleon. In terms of electric and magnetic EMFF's, 
and G^, i.e. considering the matrix element given in eq. ([T|) and the definitions of eq. ([2]), 
we have 

|GcW)l = (W)| 2 + ^|GW)| 2 ) 1/2 (l + 2 -f-)' 1/2 , 

and this is the relation that we use to fit the data on \G?L\ for both proton-antiproton and 
neutron-antineutron production. 

The proton-antiproton production experiments were of two types, 1) the exclusive pair pro- 
duction |13H20| and 2) production of the pair with a photon [21]. In the latter case the pair 
production energy is obtained by assuming that the photon was produced by the electron or 
positron and that no other photons were emitted but undetected. Fits of the model were made 
both with and without the latter data [21j . In Figs. [5lfT2l the fit curves corresponding to the two 
possibilities: with and without BaBar data, are shown as solid and dashed lines, respectively. 
The free parameters of this model are: 

• the three cut-offs: Ai, A2 and A^> which parameterize the effect of hadronic FF's and 
control the transition from non-perturbative to perturbative QCD regime in the 'jNN 
vertex; 

• five pairs of vector meson anomalous magnetic moments and photon couplings (k , g a /fa), 
with a = p, p' , co, u', 4>. 

The best values for these 13 free parameters together with the constants of this model are 
reported in Table [3j 

The fixed parameters concern well known measurable features of the intermediate vector mesons 
and dynamical quantities. Particular attention has to be paid to Aqcd- In fact we use the 
values Aqcd = 0.15 GeV in all cases but for the case=l without BaBar data, where instead: 
Aqcd = 0.10 GeV. The use of such a reduced value is motivated by the requirement of having 
no real poles in meson-nucleon and quark-nucleon FF's (Sec. [3]). As Aqcd = 0.15 GeV is closer 
to the values preferred by high energy experiments, it suggests that case=s is the more physical 
model. Another reason to prefer it on physical grounds is that the width formula of the vector 
meson decay in case=s is determined by relativistic perturbation theory. Case=l was chosen 
because it is a simpler relativistic modification of the non-relativistic Breit-Wigner form. This 
in our view is a less physical reason. 



7 Discussion 

The Lomon-Gari-Krumpelman Model [5] was developed for and fitted to space-like EMFF data. 
To enable the model to include the time-like region only the vector meson (of non-negligible 
width) propagators needed revision to appropriately represent a relativistic BW form at their 
pole in the time-like region. Two such forms are discussed above, case=l, the minimal alteration 
from the non-relativistic BW form, and case=s derived from relativistic perturbation theory. 
The resulting modification in the space-like region is minor and affected the fit there very little. 
With the new form of the vector meson propagators the simultaneous fit to the space-like EMFF 
and the time-like nucleon-pair production data was satisfactory as seen in Figs. HlfTTI and by the 
X 2 values of Table [3j 

The x 2 contributions from each space-like EMFF differ little between case=l and case=s and 
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Parameter 


case = s 
With BaBar 


case = 1 
With BaBar 


case = s 
No BaBar 


case = 1 
No BaBar 


9p/f P 


2.766 


2.410 


0.9029 


0.4181 




-1.194 


-1.084 


0.8267 


0.6885 


M p (GeV) 


0.7755 (fixed) 


T p (GeV) 


0.1491 (fixed) 


dull fu) 


-1.057 


-1.043 


-0.2308 


-0.4894 




-3.240 


-3.317 


-9.859 


-1.398 


M u (GeV) 


0.78263 (fixed) 


9<j>/f<j> 


0.1871 


0.1445 


0.0131 


0.1156 


K<f> 


-2.004 


-3.045 


37.218 


-0.2613 


M (GeV) 


1.019 (fixed) 


H (GeV) 


20.0 (fixed) 


9u'/fu> 


2.015 


1.974 


1.265 


1.649 




-2.053 


-2.010 


-2.044 


-0.6712 


M u i (GeV) 


1.425 (fixed) 


T^GeV) 


0.215 (fixed) 


9 P '/f P ' 


-3.475 


-3.274 


-0.8730 


-0.0369 




-1.657 


-1.724 


-2.832 


-104.35 


Mpi (GeV) 


1.465 (fixed) 


T P ' (GeV) 


0.400 (fixed) 


Ai (GeV) 


0.4801 


0.5000 


0.6474 


0.6446 


A 2 (GeV) 


3.0536 


3.0562 


3.0872 


3.6719 


Ad (GeV) 


0.7263 


0.7416 


0.8573 


0.8967 


Aqcd (GeV) 


0.150 


0.100 



Table 3: Best values of fit parameters and constants. 



are approximately the same as in the space- like only fit of Ref . |12| . 

However the fit in the time-like region, as measured by x 2 > is qualitatively poorer when the 
BaBar data [21] are included (x 2 /d.o.f.=2.5) than when that set of data is omitted (x 2 /d.o.f.=0.5 
for case=l, and is 1.0 for case=s). As the quality of the fit is poorer when the BaBar data 
are included it may indicate an inadequacy in the model. However the energy of the nucleon 
pairs produced in the BaBar experiment, unlike that of the exclusive pair production [13]- [20] . 
depends on the assumption that the observed photon is from electron or positron emission and 
is not accompanied by a significant amount of other radiation. The resultant theoretical error is 
not fully known although relevant calculations have been made [22]. The angular distributions 
may be sensitive to these radiation effects affecting the values of the \G E /G P M \ ratio whose data 
are displayed in Fig. [12] together with our prediction. 

Figures @][9] are extended to higher momentum-transfers than the present data to show how the 
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Figure 4: Space-like magnetic proton EMFF Figure 5: Space-like electric proton EMFF nor- 
normalized to the dipole and fi p , in case=l and malized to the dipole, in case=l and case=s, 
case=s, including and not the BaBar data. including and not the BaBar data. 
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0.1- 
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Figure 6: Space-like magnetic neutron EMFF Figure 7: Space-like electric neutron EMFF 
normalized to the dipole and fi n , in case=l and normalized to the dipole, in case=l and 
case=s, including and not the BaBar data. case=s, including and not the BaBar data. 



four different fits may be discriminated by new data. Figure [8] for R p indicates that at the higher 
momentum-transfers extended data may discriminate the smaller case=s no BaBar prediction 
from the larger case=l and case=s wMi-BaBar predictions and from the still larger case=l 
iio-BaBar prediction. Figure [9] for R n shows that at high momentum-transfer the case=s pre- 
dictions are higher than those for case=l. 
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Figure 8: Space-like ratio R p = fi p G p E /G p RI 
normalized to fi p , in case=l and case=s, in- 
cluding and not the BaBar data. 
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Fi gure 9: Space-like ratio R n — /i n G E j G 7 ^ 
normalized to fi n , in case=l and case=s, in- 
cluding and not the BaBar data. 
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Figure 10: Time-like effective proton FF data Figure 11: Time-like effective neutron FF data 
(nine sets |13H21| ) and fit, in case=l and (only FENICE [13]) and fit, in case=l and 
case=s, including and not the BaBar data. case=s, including and not the BaBar data. 



Figure [TT] is extended in energy for the same reason. It clearly shows that at higher energy 
case=l ho-BaBar may be discriminated from the other three fits by moderately precise data. 
An extension of Fig. [10] would only show the production of proton pairs remaining very close to 
zero. However in the range of energy already covered it is evident that the case=s iio-BaBar 
result is difficult to reconcile with the BaBar data for s = 5 — 7 GeV 2 . However for s = 7.5 — 8.5 



15 



Case=s with BaBar 
Case=l with BaBar 
Case=s no BaBar 
- Case=l no BaBar 




8 9 10 
q 2 (GeV 2 ) 



K3 



^5 

o 

Q.BJ 

o 



— Case=s with BaBar 

- Case=l with BaBar 

2 - - Case=s no BaBar 
Case=l no BaBar 




10 



q 2 (GeV 2 ) 



Figure 12: Modulus of the ratio G P E /G P M , Figure 13: Prediction for the phase of the ratio 

data [21] and prediction, in the time-like re- G P E /G P M , in the time-like region, in case=l and 

gion, in case=l and case=s, including and not case=s, including and not the BaBar data, 
the BaBar data. 



GeV 2 the ho-BaBar fits are closer to the BaBar data than are the wMi-BaBar fits. 
Figures PT21 and [T3l show that experiments in the time-like region for the ratio and 
the phase difference of G P E and G P M would be effective in discriminating between the models 
presented here and other models as well. 
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Appendix A: 
Dispersion Relations 

Dispersion relations are based on the Cauchy theorem. Consider a function F(z), analytic in 
the whole z complex plane with the discontinuity cut (so,oo). If that function vanishes faster 
than 1/ln \z\ as \z\ diverges we can write the spectral representation 

7T 



This is the so-called DR for the imaginary part where it is understood that the imaginary part 
is taken over the upper edge of the cut. 

The extension to the case where there is a finite number of additional isolated poles is quite 
natural. Indeed, considering a function with the set of N poles {zj} (J = 1, . . . , N) of Sec. 14. 1\ 
under the same conditions we obtain the spectral representation 



1 f°° Im[F(x)]dx 



[F(z') 1 If 
F(z)+2niJ2 Res Zj =- 

j=l L J Js t 



i so X z 



(A.2) 



where Res[<?(z), Zq\ stands for the residue of the function g(z) at z = zq. Furthermore, since we 
know the poles, we can use the more explicit form 



/(*) 



F{z) = rpv ( r 

Uj=i(z - Zj) 

where f(z) is the pole-free part of F(z), but it has the same discontinuity cut. Using this form 
in the residue definition of eq. (|A.2p and defining F(z) as the regularized version of F(z), we 
have 

tZi Uk=l,k*j( z k -zj)Zk-z vr J SQ x-z 

which is exactly the same expression as eq. (jllj) . In other words, the DR procedure, using 
only the imaginary part of a generic function, which is suffering or not from the presence of 
unwanted poles, guaranties regularized analytic continuations, the poles, even if unknown, are 
automatically subtracted. 

Appendix B: 
A third case 

We consider a regularized vector meson propagator [23] 

1 



D(s) 



Ml - s + n(s) ' 



where Mq is the bare mass of the meson and H(s) is the scalar part of the tensor correlator. The 
imaginary part, due to the pion loop, can be obtained using the so-called Cutkosky rule [23] as 



i m n W = -W^*( s - So ), 7o - fff^, . (B.i) 



17 



The real part of II(s) represents the correction to the bare mass Mq in such a way that the 
dressed mass becomes 

M 2 = M 2 + Ren(s). 

It follows that the propagator can be written in terms of M 2 and the only imaginary part of 



D(s) 



1 



1 



M 2 -s-iImn(s) M 2 - s — i 9{s - S ) 70 y/(s - s ) 3 /s ' 



(B.2) 



Actually, only the imaginary part of this expression makes sense because of the Heaviside step 
function in the definition of eq. (|B.ip . nevertheless, using DR, one can determine the complete 
propagator starting just from its imaginary part. The propagator is expected to be real below 
the threshold sq. In particular, using eq. (|A.1|) for t < sq, we have 



7T 



so 



lmD(s)ds 70 
s — t TV 



so 



\ 3 ds 



[ s (M2- S )2 +7o 2( s _ s ~ )3] (s _ t) 



while the real part over the time-like cut (so,oo), i.e. for s > sq, is 

lmD(s')ds' 



i 

ReD(s) = -Pr, 

7T 



70 Pr 

7T 



^'(s'-so) 3 ^' 



(B.3) 



(B.4) 



, So [ s '(M2- s 2 +7oV-*o) 3 ]( S '-s)' 

In this case the "natural" space-like extension of the original form given in eq. (IB.2[) is no more 
possible, in fact such a form, when we forget the Heaviside function in the denominator, develops 
a second cut which extends over the whole space-like region. It follows that we can not write 
an expression like 



3 D(s < s ) / 



A/ 2 



E 



R k 



s k 



Physical poles 



where we get, in the space- like region, a regular and real propagator simply by subtracting the 
physical poles. 

The only possibility to go below threshold is to use the DR's of eq. (|B.3p and (|B.4p . We compute 
explicitly the DR integrals using the substitution 

so 



x 



1 



ds = 2sq- 



x dx 



(1-x 2 ) 2 

S G (Sq, OO) — > X G (0, 1) 



The regularized form for D(s) is 



D(a) 



3 f 3 In ( 6±i 



4) 



s < s 



I go-i I 



i=l 



7TS70 



■ + 



iloVs(s-so) 3 
s(M 2 -s) 2 + 1 2 (s-s f 



s > s 



18 



The four values x\ (i = 0, 1, 2, 3), with = y x 2 , are the roots of the 4 th -degree polynomial in 
x 2 , which represents the denominator of the integrands in both DR's: 

{ [(1 - x 2 ) M 2 /s - l] 2 h 2 + x 6 } (a„/t - 1 + x 2 ) , (B.5) 

in particular: x 2 , = 1 — sq/s, is the only root that depends on s, while the three Xi, with i = 1, 2, 3, 
are the constant zeros of the first polynomial factor of two in eq. (11151) . The value at s = can 
be obtained as 



5(0) = 5— V 



Infill 

\Xi-l 



Concerning the asymptotic behavior, when s — > ±oo, i.e.: Xq — > 1, is 

1 In |s| 7q ln|s| 



Da) 



\s\^ TTJoUUli 1 ~ A) S 7r ( 1 +7o) S 



where the last identity follows because the product at denominator is just the 3 th degree x 2 - 
polynomial of eq. (|B.5p evaluated at x 2 = 1. 

A data fit was not made for this case because the resonance shape it produces is intermediate 
between the fitted case=l and case=s. 



Appendix C: 

The threshold behavior 

The effective proton and neutron EMFF's extracted from the cross section data through the 
formula of eq. (|17p have a quite steep enhancement towards the threshold, i.e. when q 2 — > 
(2Mtv) 2 . This is a consequence of the almost flat cross section measured in the near-threshold 
region: (2M^) 2 < q 2 < (2 GeV) 2 . Such a flat behavior is in contrast with the expectation in case 
of a smooth effective FF, which gives, near threshold, a cross section proportional to the velocity 

of the outgoing nucleon \J ^ ~ / q 2 . Moreover, in the threshold region the formula of eq. (|17p 

has to be corrected to account for NN finale state interaction. In particular, in the Born cross 
section formula, in case of proton-antiproton, we have to consider the correction due to their 
electromagnetic attractive interaction [25 j . Such a correction, having a very weak dependence 
on the fermion pair total spin, factorizes and, in case of pointlike fermions, corresponds to the 
squared value of the Coulomb scattering wave function at the origin, it is also called Sommerfeld- 
Schwinger-Sakharov rescattering formula [26] . Besides the Coulomb force also strong interaction 
could be considered. Indeed, when final hadrons are produced almost at rest they interact 
strongly with each other before getting outside the range of their mutual forces |27| . Indeed 
there is evidence for near threshold quasi-bound iViV states with widths in the tens of MeV [28] . 
It follows that EMFF values in this energy region are affected by different kinds of corrections 
whose form and interplay are not well known. Hence we decided to include in the present 
analysis only data above q 2 = 4 GeV 2 , to avoid the threshold region. 

Figures [2] and [15] show the residue data-over-fit for the proton and neutron effective FF's, 
respectively. They have been obtained dividing the fit functions shown in Figs. [10] and [TT] by 
the corresponding data on |G^|. The threshold enhancement of the proton data exceeds the 
fit by a factor of more then two and, in the neutron case, even within large errors, the factor is 
about three. 
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Figure 14: Residue for the proton effective FF: Figure 15: Residue for the neutron effective FF: 

I G cS, data I / 1 G cS, fit I . where I G cff, fit I haS been ob " I G eS, data I / 1 G cff, fit I ' where I G eS, fit I haS been ob- 
tained considering only data with q 2 > 4 GeV 2 . tained considering only data with q 2 > 4 GeV 2 . 



References 

[1] C. F. Perdrisat, V. Punjabi and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 59, 694 (2007) 
|arXiv:hep-ph/0612014] . 

[2] F. Iachello, A. D. Jackson and A. Lande, Phys. Lett. B 43, 191 (1973). 

[3] M. F. Gari and W. Kruempelmann, Phys. Lett. B 274, 159 (1992) [Erratum-ibid. B 282, 
483 (1992)]. 

[4] See e.g. E. Tomasi-Gustafsson, F. Lacroix, C. Duterte and G. I. Gakh, Eur. Phys. J. A 24, 
419 (2005) iarXiv:nucl-th/0503001|. 



[5] E. L. Lomon, |arXiv:nucl-th/0609020t 

E. L. Lomon, Phys. Rev. C 66, 045501 (2002) [arXiv:nucl-th/0203081| ; 

E. L. Lomon, Prepared for 9th International Conference on the Structure of Baryons 

(Baryons 2002), Newport News, Virginia, 3-8 March 2002; 

E. L. Lomon, Phys. Rev. C 64, 035204 (2001) [arXiv:nucl-th/0104039| . 

[6] L. L. Foldy, Phys. Rev. 87, 688 (1952). 

[7] L. N. Hand, D. G. Miller and R. Wilson, Rev. Mod. Phys. 35, 335 (1963). 

[8] P. Merge!!, U. G. M eissner and D. Drechsel, Nucl. Phys. A 596, 367 (1996) 
[arXiv:hep-ph/9506375] . 

[9] S. J. Brodsky and G. R. Farrar, Phys. Rev. Lett. 31, 1153 (1973); 
V. Matveev et al, Nuovo Cimento Lett. 7 (1973) 719; 
G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 



20 



[10] R. J. Eden, P. V. Landshoff, D. I. Olive, and J. C. Polkinghorne, The Analytic S-matrix 
(Cambridge University Press, Cambridge, England, 1966); 
G. F. Chew, The Analytic S-Matrix (Benjamin, New York, 1966). 

[11] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[12] C. Crawford et al, Phys. Rev. C 82, 045211 (2010) [arXiv: 1003.09031 [nucl-th]]. 

[13] A. Antonelli et al. [FENICE Collaboration], Nucl. Phys. B 517, 3 (1998). 

[14] B. Delcourt et al. [DM1 Collaboration], Phys. Lett. B 86, 395 (1979). 

[15] D. Bisello et al. [DM2 Collaboration], Nucl. Phys. B 224, 379 (1983). 

[16] M. Ablikim et al. [BES Collaboration], Phys. Lett. B 630, 14 (2005) [arXiv:hep-ex/0506059] . 

[17] T. K. Pedlar et a l. [CLEO Collaboration], Phys. Rev. Lett. 95, 261803 (2005) 
[arXiv:hep-ex/0510005] . 

[18] G. Bardin et al. [LEAR Collaboration], Nucl. Phys. B 411 (1994) 3. 

[19] T. A. Armstrong et al. [E760 Collaboration], Phys. Rev. Lett. 70, 1212 (1993). 

[20] M. Ambrogiani et al. [E835 Collaboration], Phys. Rev. D 60, 032002 (1999). 

[21] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 73, 012005 (2006) 
[arXiv:hep-ex/0512023] . 

[22] V. V. Bytev, E. A. Kuraev, E. Tomasi-Gustafsson and S. Pacetti, Phys. Rev. D 84, 017301 
(2011) [arXiv: 1103.44701 [hep-ph]]. 

[23] F. Klingl, N. Kaiser and W. Weise, Z. Phys. A 356, 193 (1996) [arXiv:hep-ph/9607431] , 

[24] C. Itzykson and J.-B. Zuber, Quantum Field Theory, Mc Graw-Hill, New York (1980). 

[25] R. Baldini, S. Pacetti, A. Zallo and A. Zichichi, Eur. Phys. J. A 39, 315 (2009) 
|arXiv:071 1.17251 [hep-ph]]. 

[26] A. D Sakharov, Zh. Eksp. Teor. Fiz. 18, 631 (1948) [Sov. Phys. Usp. 34, 375 (1991)]; 

A. Sommerfeld, Atombau und Spektralliniem (Vieweg, Braunschweig, 1944), Vol. 2, p. 130; 
J. Schwinger, Particles, Sources, and Fields, Vol. Ill, p. 80. 

[27] J. Haidenbauer, H. W. Hammer, U. G. Meissner and A. Sibirtsev, Phys. Lett. B 643, 29 
(2006) [arXiv:hep-ph/0606064] . 

[28] S. Wycech and B. Loiseau, arXiv-hep-ph/0508064vl(2005). 



21 



